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Abstract 

Topologically non-trivial superconductivity has been predicted to occur in superconductors with 
a sizable spin-orbit coupling in the presence of an external Zeeman splitting. Two such systems have 
been proposed: (a) ,s-wave superconductor pair potential is proximity induced on a semiconductor, 
and (b) pair potential naturally arises from an intrinsic s-wave pairing interaction. As is now well 
known, such systems in the form of a 2D film or ID nano-wires in a wire-network can be used 
for topological quantum computation. When the external Zeeman splitting T crosses a critical 
value Fc, the system passes from a regular superconducting phase to a non-Abelian topological 
superconducting phase. In both cases (a) and (b) we consider in this paper the pair potential 
A is strictly s-wave in both the ordinary and the topological superconducting phases, which are 
separated by a topological quantum critical point at Tc = ^/^~+~|J?, where n{» A) is the 
chemical potential. On the other hand, since Fc >> A, the Zeeman splitting required for the 
topological phase (F > Fc) far exceeds the value (F ~ A) above which an s-wave pair potential is 
expected to vanish (and the system to become non-superconducting) in the absence of spin-orbit 
coupling. We are thus led to a situation that the topological superconducting phase appears to 
set in a parameter regime at which the system actually is non-superconducting in the absence of 
spin-orbit coupling. In this paper we address the question of how a pure s-wave pair potential can 
survive a strong Zeeman field to give rise to a topological superconducting phase. We show that 
the spin-orbit coupling is the crucial parameter for the quantum transition into and the robustness 
of the topologically non-trivial superconducting phase realized for F >> A. 

PACS numbers: 03.67.Lx, 71.10.Pm, 74.45. +c 
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I. INTRODUCTION 



Recently topologically non-trivial superconductivity has been theoretically predicted to 
occur in two classes of systems with spin-orbit (SO) coupling. They are: (a) SO coupled 
semiconductors in which s-wave superconducting pair potential is induced by the proximity 
effectP^ and (b) SO coupled systems with superconductivity due to intrinsic s-wave pairing 
interactionSEl. A third system - surface of a 3D strong topological insulator (TI) - can also 
support topological superconductivity when the latter is proximity induced^. In this paper 
we will ignore the latter and concentrate only on (a) and (b). In both classes (a) and (b) 
the SO coupling is a consequence of the breakdown of the structural space inversion (SI) 
symmetry. We will take the resultant SO coupling to be of the Rashba type. A system 
in class (a) can be artificially grown as a heterostructure consisting of a SO coupled semi- 
conductor in proximity contact with a s-wave superconductor^. An example of class (b) 
is a non-centrosymmetric superconductor'^ or an s-wave Feshbach resonant system with an 
accompanying spin-orbit coupling and Zeeman splitting both of which can be created in a 
cold fermion atomic system^ ^. There is also another class (class (c)) of topological super- 
conductors with non-Abelian statistics, which have been studied extensively in the recent 
literature. Here the topological nature arises entirely from intrinsic chiral p-wave super- 
conductivity without having any underlying s-wave superconductivity in the system. Some 
examples of this class (c) are the even-denominator 5/2 fractional quantum Hall effect!^, 
superconducting strontium ruthenat^, A-phase of superfluid He-Jp^l^ and superfluid ultra- 
cold fermionic gases based directly on the p-wave Feshbach resonanc^^. These class (c) 
non-Abelain superconductors are effectively equivalent to being spinless, i.e. completely 
spin-polarized, and are therefore immune to any Zeeman splitting to the leading order. We 
do not discuss the class (c) systems in this paper since the main conceptual issue being 
addressed in this paper does not apply to these systems. 

The topologically non-trivial superconducting systems mentioned above are characterized 
by order parameter defects, such as a vortex and a sample edge, which carry a unique 
bound state at zero excitation energjff^. These bound zero energy modes, called Majorana 
fermion modes after E. Majorana?^, can actually be thought of as particles which are their 
own anti-particle^. In other words, they are represented by second quantized operators 7 
which satisfy the hermiticity condition 7^ = 7. This is strikingly different from the regular 
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FIG. 1. (Color online) The two spin-orbit bands with (blue curves) and without (red curves) 
Zeeman splitting are shown schematically. With Zeeman splitting the bands have a band gap at 
the origin. When the Zeeman splitting is large enough so that the chemical potential (dotted circle) 
lies in the gap, the system has only one Fermi surface and the ordinary superconducting phase 
gives way to a topologically non-trivial superconducting phase. 

fermionic modes for which the second quantized operators 7^ c. Because of the existence 
of the zero-energy Majorana modes, a 2D non-Abelian topological superconductor with 2N 
vortices each carrying a single Majorana mode constitutes a system with a 2^-fold ground 
state degeneracy protected by an excitation gapP^. When a Majorana fermion mode in a 
given order parameter defect is adiabatically moved (braided) around another Majorana 
mode, the initial state transforms into another one which is a different linear combination 
in the same ground state manifold. These unitary transformations within the ground state 
manifold are manifestations of the non-Abelian statistics of the Majorana fermions. The non- 
Abelian statistics of the Majorana fermions has recently come under intense focus because of 
its potential application in fault tolerant topological quantum computation (TQC)'^^'^. Both 
2D spin-orbit coupled superconducting films and its ID version as quantum nanowires in a 
wire-network have been proposed as potential platforms for TQCP^^^^. In this paper we will 
not discuss the Majorana fermion modes or their potential application to TQC. Instead we 
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will focus on the bulk superconducting phases and the quantum phase transitions (QPT) by 
exploring the interplay of the spin-orbit coupling, Zeeman splitting, and the s-wave pairing 
interactions. 

Both systems in classes (a) and (b) are in an ordinary (non-topological) superconducting 
phase in the absence of an externally imposed Zeeman splitting. The Zeeman splitting, 
which can be applied either by a parallel magnetic fielcP^ or by the proximity effect of a 
nearby magnetic insulatoi'^, creates a gap between the two spin-orbit bands as shown in 
Fig. [!} From mean field calculation^^ it is clear that when such a Zeeman splitting F crosses 
a critical value, Fc = a/A^ + /i^, where A is the s-wave superconducting pair potential and 
fi is the chemical potential, the system makes a transition to a topological non-Abelian 
superconducting phase. The critical Zeeman splitting corresponds to the value at which the 
underlying Fermi surface in the absence of superconductivity shifts from being at both SO 
bands to occupying only the lower SO band (Fig. [T]). Thus, at this value of the Zeeman 
splitting, the underlying system changes from having two Fermi surfaces (one in each band) 
to having just one. 

Since /i usually far exceeds A, the critical value of the Zeeman splitting Fc also far exceeds 
A. Therefore, in the absence of SO coupling, Fc far exceeds the Zeeman splitting (F ~ A) 
above which an s-wave pair potential should decay to zero. The loss of an s-wave pair 
potential due to a strong Zeeman sphtting is due to the fact that for F ^ A forming a spin- 
singlet pair potential with zero net momentum is impossible (for the discussion of topological 
phase transition induced by Zeeman splitting we will ignore the states with non-zero values 
of the Cooper pair momentum). Now let us emphasize the fact that in both cases we consider 
in this paper (superconducting pair potential proximity induced, and superconducting pair 
potential due to an intrinsic on-site pairing interaction), the pair potential is strictly s-wave 
in both the ordinary and the topologically non-trivial superconducting phases. The fact 
that the pair potential is s-wave when it is proximity induced from a s-wave superconductor 
is self-evident. That it remains s-wave (and is not a mixture of s- and p-wave due to the SO 
coupling) even when the pair potential is due to an intrinsic on-site pairing interaction is not 
so obvious. In this case, the pure s-wave symmetry of the pair potential follows from the fact 
that the intrinsic pairing interaction we consider is spatially local, and thus forming a p-wave 
component of the pair potential is forbidden by the fermion anticommutation relation (for 
a more detailed discussion see Sec. V). Since the pair potential is purely s-wave in both 
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classes (a) and (b) and in both phases (ordinary and topological) in each, how an s-wave pair 
potential survives a strong Zeeman splitting F >> A to realize the topologically non-trivial 
phase is the central conceptual question we address in this paper. 

The basic conceptual issue being discussed here is the topic often alluded to as the 
Chandrasekhar-Clogston (CC) limilP^'^Sl ordinary s-wave superconductivity, which states, 
in effect, that an s-wave superconductor, where the Cooper pairing is between spin- up and 
spin-down electrons near the Fermi surface, cannot withstand a Zeeman splitting larger than 
the superconducting gap. This is because then spin splitting exceeds the superconducting 
gap energy, making it impossible for a superconducting ground state to develop. On first 
sight, it appears that the condition on the Zeeman splitting needed for superconductivity 
in Refs. [1-4] far exceeds this limit, thus destroying all superconductivity! This has caused 
some confusion about the very existence of the topological superconducting phase using 
either a heterostructurd^ where s-wave superconductivity is induced by proximity effect or 
using SO-coupled systems with intrinsic s-wave pairing interaction^'. 

The mean field calculations of Ref. [1-4] are not enough to resolve this question. This 
is because a mean field theory is not just a postulate to assume the existence of a mean 
field pair potential in the Hamiltonian H as is done in these works; one is also required to 
establish the finiteness of the pair potential by satisfying the self-consistent gap equation. 
In other words, we need to satisfy the Bardeen-Cooper-Schrieffer (BCS) self-consistent gap 
equation with a strong Zeeman potential (F > a//!^ + A^) to check if a non-zero s-wave pair 
potential A gives a consistent solution. This will ensure that the mean field H in Refs. [1-4] 
is not flawed to begin with, and our BdG solution of the Majorana fermion is not a spurious 
mathematical result with no physical connection. In this paper we perform this study by 
self-consistently solving the appropriate BCS gap equations in the presence of attractive 
s-wave pairing interaction, SO coupling, and an externally applied Zeeman splitting. Note 
that satisfying such a gap equation is a requirement for superconducting pair potential when 
it is derived from the microscopic pairing interactions. However, when the pair potential 
is proximity induced on a spin-orbit coupled system by a nearby s-wave superconductor, the 
gap equation need not be satisfied. In this case, the SO coupled system simply 'inherits' the 
pair potential of the nearby superconductor. 

For the case of intrinsic pairing interactions, we show that the s-wave pair potential 
indeed remains non-zero even beyond the Zeeman splitting above which it would be lost in 
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the absence of SO coupling. In fact, in the non-Abehan phase the non-zero value of the 
pair potential crucially depends on and increases with the SO coupling constant, which is 
consistent with the fact that it is zero in the absence of the SO coupling. A simple intuitive 
way to understand this starts by recalling how s-wave superconductivity is destroyed by 
a Zeeman splitting. In the absence of SO coupling, the two spin bands are shifted by an 
energy gap proportional to the Zeeman splitting F. With increasing F it becomes increasingly 
difficult for the system to create s-wave spin-singlet pairs with zero net momentum. Finally, 
when F crosses a value ~ A the s-wave pair potential vanishes. The critical Zeeman splitting, 
Fc = y^A^~+/?, needed for the topological phase transition is thus squarely beyond the 
acceptable Zeeman splitting the pair potential can sustain. It is now important to realize 
that, in the presence of the SO coupling, the two spin-orbit bands cannot simply be viewed as 
'spin-up' and 'spin-down' bands. Instead, they both have a non-zero minority spin amplitude 
coexisting with the majority spin component. Therefore, even when the Zeeman splitting 
is large enough to make the Fermi surface lie only in the lower band, spin-singlet s-wave 
pairs cannot be completely lost. If the superconductivity is due to an intrinsic pairing 
interaction the gap equation shows that the pair potential, though always non-zero, decays 
with increasing Zeeman couplings. However, the magnitude of the pair potential in the 
non-Abelian phase can be increased by increasing the magnitude of the spin-orbit coupling, 
which, therefore, enables a stable non-Abelian phase in the phase diagram . On the other 
hand, if the superconductivity is due to the proximity effect, there is no requirement of 
satisfying the self-consistent gap equation. In this case, which applies to the heterostructure 
geometry, the superconducting pair potential is simply 'inherited' from the adjacent s-wave 
superconductor (Sec. VIII). 



II. HAMILTONIAN 

We assume that the quasi-2D electron system is described by the model Hamiltonian 

H^Ho + Hso + Hr + Hint, (1) 

where Hq describes the bulk conduction electrons, Hso is the spin-orbit interaction term, 
Hy represents the Zeeman coupling, and i^mt represents the electron-electron interaction. 
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Explicitly, we have 

p,a 

HsO = a ^ 4^ {pyT^c - PxTy)^^^, Cpa', (3) 
p,a,(r' 

= r J] (^cj^cpt - cJ,^Cp^) , (4) 
p 

^int = ^ X] Z]^(^)^I'+9<^4'-9^'V^'Cp,<T, (5) 
p,p',q a,a' 

where — p^/2m—^ is the bulk spectrum (measured relative to the chemical potential //), a 
is the strength of the Rashba spin-orbit coupling, F represents the Zeeman field, V{q) is the 
short-ranged interaction potential (we will later restrict ourselves only to an on-site pairing 
interaction, V{q) independent of q, natural for s-wave order), Tx{y) are Pauli matrices, and 
Cpo- (cpcr) is the creation (annihilation) operator corresponding to the single-particle state 
with momentum p and spin a. It is convenient to work in the spinor basis provided by the 
eigenfunctions 0a (p) of the single-particle Hamiltonian Hq -\- Hso, 

Up) = 4 I V, I ' (6) 

where A = ± and e*^'' = {px + ipy) jp- The corresponding eigenvalues are 

eA(p) = ^p + Aap. (7) 

The electron c-operators can be expressed in terms of the annihilation operators ap\ associ- 
ated with the spinor eigenstates as 

^v<y = ^ <^a(p, cr)apx. (8) 

A 

Using the spinor representation, the Hamiltonian becomes 

Ho + Hso = ^M^lx^^px, (9) 

p,A 

Hr^rJ^^U^P-x (10) 



p,A 

2 



Hint J2 ^(?)XAA'(P + 9, p)Xm/.'(p' - Q, P'Wp+qxO,l,'-g^,ap'^,'apx, (11) 



P,P',9 A,A',/i,//' 

where 



XAA'(Pi,P2) = ^ (1 + AA'e-(V-^-)) (12) 



is the scalar product of two spinors 0a (Pi) ^-nd 0a'(P2)- 



III. SUPERCONDUCTING GAP EQUATIONS 



To derive the gap equations, we first introduce the regular and anomalous Green functions 



as 



Qxx'{p,T- t') = -{Trapx{T)aly{T')), 
J^xx'iP,r- t') = \{Trapxir)a^pX'ir')), 
7tAp.r- t') = \{Tral,iT)al^y{T')), 



(13) 
(14) 
(15) 



where is the time ordering operator and the operators apx{T) are in the Heisenberg 
representation. The correlation functions J^xx'ip) = J^xx'{p,0+) have the properties 



J^xx{-p) = -J'aa(p), 
J^+-i-p) = J^-+ip), 



J-+(-p)= J-,\(p), 



(16) 
(17) 



The definitions (14) and (15) of the anomalous correlation functions follow the convention 



r — r 



used by Gor'kov and RashbaP^. 

Following the standard procedure, we write the equations of motion for the Green func- 
tions using the time evolution of the a-operators, drapx^r) = [H,apx]. We have 

[-•9, - exip)]gxx'{p,r- r') - rg_xx'ip,r - r') 

+ ^(P-9)XAM(P,g)XAi/.i(-p, -g)Ai/iiJ'^i^.(-9>0+)-^AiA'(-P>^-7-') =hx'S{i 

q,fj.,fj.l,Xi 

[-dr + exipWU-P, T-r')- m,,,(-p, r - r') (18) 
- A XI "^^P ^ g)XMA(-g, -p)XmAi (g, P)/^'^wi (-9, 0+)+^Aia(-P, t-t')= 0. 

qr,Ai,/^l,Ai 

The gap function can be defined as, 

AaAi(p) = Ai X ^(P- g)XA^XP>g)XAiMi(-P>-9)/^i-^mM(-9'0+)- (19) 



Introducing the definition of the gap function in Eq. (18), we have 



r — r 



[-dr - ex{p)]Gxx'{p. T-r')- TG_xx'{p. r - r') + ^xxAp)Kx'i-P^ r - r') = 5,v5(- 
[-dr + 6,(p)]J-+,(-p,r - r') - r^+,,,(-p,r - r') + > ^ Al,,^,,,(p, r - r') = 0. (20) 



m,,,(-p,r - r') + ^ Al^,^,,,(p,r - r') = 0. 

Ai 



Defining the Fourier transforms of the correlation functions in the usual way, Qxy {p, t) = 
ksT e~'"^"'^Qx\'{p, i^n)^ the set of equations of motion can be expressed in a matrix form 
as 



(21) 



where the arguments of the Green functions have been omitted for simplicity. A similar set 
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of equations, which can be obtained from (21 ) by switching the + and — labels, couples Q , 



T , and T^-. The superconducting spectrum can be obtained from the condition that 



the determinant of the 4x4 matrix in Eq. (21 ) vanish taking iiOn — )■ E. Also, by solving the 



system of equations of motion for and introducing the solutions in Eq. ( 19 ) we obtain 
the self-consistent gap equations. In general, we have Aaa(p) = «e~*^p[Aos(p) + AAoa(p)], 



where the symmetric and antisymmetric components of the diagonal gap functions are 
Ao.(p) = 



Aoa(p) = -iY. + [( - cos(^^p - Q^)e''- + z( - sin(^^p 



Similarly, the off-diagonal gap functions can be expressed as A^ -a(p) = "^e *^p[Ais(p) + 
AAi„(p)], with 



Ala(p) 



y{p-q)-y{p'rq) 



- sin(^p - eq)e*^'' + - cos(ep 



Note that Ajs(— p) = Ajs(p) and Aja(— p) = — Aja(p), i.e., Aj^ and Aja represent the 
singlet and triplet components of the gap functions, respectively. 

IV. ASSUMPTION OF LOCAL INTERACTION 



Instead of solving the complicated coupled set of gap equations above, we simplify mat- 
ters by considering the case of strictly local interactions. In other words, we neglect the 
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momentum dependence of the interaction potential, V{p) = Vq < 0. Then the only non- 
vanishing component of the superconducting gap is the singlet component, Aq^ = A, and 
it becomes momentum-independent. Since by a Zeeman splitting the singlet component of 
the gap function will be the most affected, we can make this approximation to examine the 
fate of the superconducting condensate with increasing Zeeman potential. 

For a strictly local attractive interaction, the superconducting spectrum is given by. 



(26) 



where ak = ak. Solving the kinetic equations for and J-' and using Eq (22) we obtain 

the gap equation for the strictly local attractive interaction. 



+ 



1 



+ Ef ul + El /q^2 + r2(q + |A|2) V^n + Ei ul + E, 



1 



(27) 



Taking the zero temperature limit and performing the summation over the frequencies we 
obtain 



E 



1 



+ 



1 



2E,{q) 2E2{q) /q^2 + r2(q + | A|2) \'2E,{q) 2^2(9 



1 



1 



(28) 



V. ANOMALOUS CORRELATION FUNCTIONS AND GAP FUNCTIONS IN 
THE C OPERATOR REPRESENTATION 

To obtain a deeper understanding of the singlet-triplet mixing in superconductors with 
spin-orbit couplin^^, it is useful to determine the expressions of the anomalous correlation 
functions and of the gap functions in terms of the original electron operators. We first express 
the c-operators in terms of a-operators, Cpi- = (ap+ + apJ)/\/2 and Cpj, = — ie*^p(ap+ — 
ap_)/v^, and we obtain for the singlet and triplet anomalous correlation functions the 
expressions 

{T,c.p^{r)cp^{T')) = '-e''- [-F++(-p, r - r') + J'-i-P, r - r') + F+-{-p, r - r') + -F_+(-p, r - r')] , 
{T,c.p,{r)cp,{T')) = ie(i--)^^- [-^++(-p, r - r') - J^^^{-p, r - r') + aT+_{-p, r - r') - aJ^_+{-p, r - 
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where the T anomalous functions are given by equations (14) and (15). In the hmit of local 
interactions we can determine the explicit dependence of the T functions on the parameters 
of the model using Eq. (21), and we have {J-\\i + J^_a-a') oc Ao5(p) and — J--\-\i) oc 
iape^''^pAos(p). Consequently, in the c-operator representation both the singlet and the 
triplet components of the anomalous correlation function are proportional to the s-wave gap, 
{TrC-p^{T)cp^{T')) oc Aos(p) and {TrC-p^{T)cpa{T')) OC iQ;Aos(p)(Px -icrpy), respectively. 

We emphasize that, in the limit of local pairing interaction, the anomalous correlation 
function in the c-operator representation has both singlet and triplet components, but the 
corresponding gap function is purely s-wave. To show this property explicitly, we can re- 
derive the gap equations in the c-operator representation and, instead of Eq. (20), we 
obtain 

(-5r-ep-(7r)C,(p,T-T')-(TapG'_,<,.(p,T-T')+J]A,,,(p)(T,cJ,^,ct_p,^) =5,,.5(t-t'), 

(30) 

where Gaa'{p,T — t') — ~{TTCpa{T)Cp^,{T')) is the normal Green function and the gap 
functions are defined as 

Ka'ip) ^-J2^iP- q){c-pa{r)cpAr)). (31) 

The equal time anomalous correlation {c^pa-{T)cpa'{T)) can be expressed in terms of Txy{0). 
Explicitly, we have 

9 9 

In the limit of strictly local interactions V{p—q) — V{p+q) — Vq and oc qx—iqy 

and, consequently, the triplet component vanishes, A-|-|- = = 0. On the other hand, the 
expression for the singlet component of the gap becomes identical with the right hand side 
of Eq. (22), hence we have A^^{p) = Aos(p) = A. We conclude that in a superconductor 
with spin-orbit coupling and on-site pairing interactions the anomalous correlation function 
is characterized by a mixture of singlet and triplet components, yet the gap function has 
purely s-wave symmetry. A p-wave component of the gap can develop only in the presence 
of non-local pairing interactions. As a consequence, in a system with strictly local pairing 
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interaction, if the singlet anomalous correlation vanishes, the superconducting gap as well 
as all the other components of the anomalous correlation function will vanish. 



VI. ANALYSIS OF THE GAP EQUATION 



We first analyze the gap equation, Eq. (28), in some special cases for which the solutions 



are well known. This will serve as a test for the validity of our analytical calculations. By 
putting r = and a = 0, which corresponds to the standard BCS case of a local attractive 
interaction with no spin-orbit coupling and Zeeman splitting, we find Ei = E2 = \ Q + [A^. 



In this case, from Eq. (28) we recover the standard BCS gap equation. 



E 



1 



lAI 



(32) 



where the summation over q should be performed over states satisfying < ud, with ud 
some cut-off Debye energy scale. As is well knowrpH^ since the integral on the right hand 
side (R.H.S.) diverges in the limit |A| 0, in this case a non-zero solution for A exists for 
any Vq < 0. 



Next we take the system with F = 0, a 7^ 0. In this case -Ei(2) = V [iq T- oiqY + |A| 



The gap equation now becomes. 



4 



+ |A| 



a,)2 + |A| 



(33) 



As in the previous case, the integrals on the R.H.S. of Eq. (33) diverge when |A| — )• 0, hence 
a non-vanishing solution for A exists for any Vq < 0. 

To establish the familiar result that s-wave superconductivity is destroyed by Zeeman 
splitting (in the absence of spin-orbit coupling), we consider the special case, F 7^ 0, and 
a = 0. In this case £'1(2) = 



-|- I Ap =1= F|, and the gap equation becomes 



-V. 



1 



e + w 



(34) 



where the summation over q is done over states satisfying \^q\ < ujd and + |Ap > F 



The second constraint results from the cancelations of two terms from Eq. (28) that diverge 



in the limit | A| — )• 0. Since the integral on the R.H.S. of the gap equation no longer diverges, 
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a non-zero solution for A exists only for |Vo| larger that a critical value. This implies that 
for a given strength of the attractive potential |Vo|, no non-zero solution for A can be found 
above a critical value of the Zeeman potential F. 

Finally we consider the most general case of a non-zero Zeeman potential as well as a 



non-zero spin-orbit coupling, F 7^ 0, a 7^ 0. The gap equation is then given by Eq. (28). The 
exact cancelation of the divergent terms that characterizes the a = case does no longer 



hold and the R.H.S. of Eq. (28 ) becomes arbitrarily large in the limit | A| — t- 0. Consequently, 
a non- vanishing solution for A exists for any negative value of Vq. This implies that A does 
not vanish for any value of F, or, in other words, the pair potential cannot be completely 
destroyed by a Zeeman splitting in the presence of a non- zero spin-orbit coupling. This is in 
agreement with a similar result derived earlier in a different contextP. Nonetheless, at large 
values of F the superconducting pair potential decreases exponentially with the strength of 
the Zeeman splitting as we show in the next section. 



VII. NUMERICAL SOLUTION AND QUANTUM PHASE TRANSITIONS 



Next, we determine the general solution of the gap equation by solving Eq. (28) nu 



merically. We address two distinct cases: 1) High carrier concentration regime, when the 
chemical potential fi (i.e., the Fermi energy in our zero temperature limit) represents the 
largest energy scale in the problem, /i 3> Wd, F, a/ci?, A(0), and 2) Low carrier concentra- 
tion, when ud > fi, T, akp, A{0). Here ud is the analog of the Debye frequency (i.e., the 
characteristic energy cut-off for the intrinsic pairing interaction), akp is the strength of the 
spin-orbit interaction at the Fermi wave vector, and A(0) is the value of the superconduct- 
ing gap at zero Zeeman splitting. The zero field gap is a measure of the pairing interaction 
strength, and in fact one can use Vq, instead of A(0), as an independent parameter. The 



Debye frequency acts as a cutoff in Eq. (28), i.e., the summation over g of a function 



f{q, Ei[q)) is restricted to the values of the wave vector satisfying Ei{q) < ud- 

To obtain the self-consistent numerical solution for the gap, we define the function 



1 1 

+ 



F^ 



2E^{q) 2E2{q) /q^2 + r2(q + | A|2) \'^E,{q) 2^2(9) 



1 



1. 



(35) 



With this notation, Eq. (28) becomes ^(A) = 0. This equation is characterized by two 
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FIG. 2. (Color online) The behavior of the function ^(A) (Eq. (35)) with A. The values of A 
at which ^(A) = are the solutions of the self-consistent gap equation, (a): 9{A) with A for 
various values of the Zeeman splitting T for a small value of the spin-orbit coupling constant a 
{2akF = 0.2meV) and large fj, = 0.5eV. Initially for smaller values of the Zeeman splitting there 
is just one solution for A. With increasing values of F, there are three solutions for A, one of 
which (the one in the middle) is unstable. The lower (higher) solution for A is the high (low) field 
solution which co-exist in a region of co-existence. With F going up, the high field solution becomes 
smaller but is never zero. Consequently, there is no true first order transition in the presence of a 
small a, even though there is a first-order-like precipitous drop of A at some value of F, which is a 
remnant of the first order transition for a = 0. (b): 0(A) with A for increasing values of F with a 
fixed large /i = 0.5eV and a larger spin orbit coupling lakp = OAmeV. 0(A) is now monotonically 
decreasing and the resultant unique solution for A decreases for increasing values of F. The rate 
of decrease of this solution with F is much slower and continuous. As with the case of smaller a, 
there is no quantum phase transition. 

qualitatively different regimes that are controlled by the relative strength of the spin-orbit 
interaction and the zero-field gap. If 2akF ^ ^(0), ^(A) is a monotonically decreasing 



function that starts from large positive values at A — and Eq. (35) has always one 



non- vanishing solution. By contrast, when 2akF < A(0), the function ^(A) becomes non- 
15 



monotonic for certain values of the Zeeman field F, which means that Eq. (35) can have 
multiple non vanishing solutions for a given set of parameters. To illustrate this situation, we 
show in Fig. [2] the function ^(A) for a system with large carrier concentration (/i = 0.5eV) 
and extremely low spin-orbit interaction {a = 0.5meV-Ain panel (a) and a = ImeV^-Ain 
panel (b), i.e., 2akF = 0.2meV and 2akF = OAmeV, respectively). The Debye frequency 
is ojd = 25meV and the zero field gap is A(0) = 0.5meV. At low Zeeman splitting, 
^(A) vanishes at a single point A ^ 0.5meV, but increasing F leads to a local minimum 
in ^(A) that goes to zero for F ^ 0.318meV^. Further increasing the Zeeman splitting 
leads to three non- vanishing solutions Ai(F) > A2(F) > A3(F) (see Fig. |2]), where Ai(F) 
and A3(F) are the "low field" and "high field" solutions, respectively, and A2(F) is an 
unstable solution. The "low-" and "high-field" solutions coexist in some range of Zeeman 
field strengths suggesting that the system undergoes a precipitous drop in A akin to a 
field-tuned first order phase transition. The coexistence region shrinks as the strength of 
the spin-orbit coupling increases (see Fig. ^p) and vanishes at a value ~ l.lmeV-A. 
We note that in real systems such as non-centrosymmetric superconductors the strength 
of the spin-orbit coupling is usually larger than this critical value and, consequently, the 
first-order-like precipitous drop in A may not be observable. In cold fermion systems the 
spin-orbit coupling constant can be used as a tuning parameter to interpolate between these 
two behaviors. 

The dependence of the solution of the gap equations on F is shown in Fig. |3j The 
coexistence region can be easily seen for a = 0.5mel^-A(green line in Fig[3|, corresponding 
to the (stable) solutions of the equation ^(A) = for the 9 function shown in Fig. [2^. 
Note the exponential decay of the "high field" solution A3 with increasing F. Practically 
the superconducting gap is neghgible (A < IfteV) for F > OAlmeV. The coexistence region 
shrinks as we approach the critical spin-orbit coupling (red line in Fig. [sjand Fig. ^p), then, 
for a > Oc the gap equation has a continuous solution A(F) that decreases monotonically 
with the Zeeman field (black line in Fig|3]). Note that at high fields, F > 2akF, the gap 
decreases exponentially. However, the energy scale for the spin-orbit coupling, 2akF, can be 
significant in realistic systems (tens of meV) and the high field regime may not be attainable, 
i.e., the gap will not vanish for any realistic value of the Zeeman field. 

The existence of a first-order-like drop in A that ends at a critical value of the SO coupling 
ttc is generic, i.e., this feature is present at any value of the carrier density. However, to 
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FIG. 3. (Color online) The solution of the gap equation (see Fig. [2]) plotted against the Zeeman 
splitting F for various values of the spin-orbit coupling constant a (in meV-K). For small values of 
a, the superconducting gap falls discontinuously with F but it is never strictly zero in the presence 
of the SO coupling. Consequently, there is only a first-order-like crossover which is a remnant of the 
true first order phase transition with F for a = 0. For larger values of a (black curve), the decay 
of the superconducting pair potential with the Zeeman splitting is much slower and continuous. It 
falls exponentially (but is never strictly zero) only for higher fields F > 2akF- 

realize a topologically nontrivial non-Abelian regime, it is necessary to satisfy the condition 
> /i^ + A^. Consequently, we study the solutions of the gap equation in the low density 
regime, when the chemical potential, the Zeeman field, the spin-orbit interaction, and the 
superconducting order parameter are comparable. In particular we address the following 
question: is it possible to realize the condition for the existence of a topologically nontrivial 
non-Abelian phase while maintaining a reasonable superconducting gap? Before presenting 
the results, we note that in the low-density regime the zero-field gap has a strong dependence 
on the chemical potential. More precisely, for a given set of parameters Vq, cod, and a, the 
zero-field gap A(0) decreases with /i. In our calculations we fix Vq at a value that corresponds 
to A(0) = OAmeV at /i = 2meV and, at lower carrier densities (i.e., lower values of /x), 
we calculate the zero-field gap using the gap equation. Also, we note that, as we vary the 
Zeeman splitting F, the chemical potential of a system with fixed carrier density n remains 
constant as long as the high-energy band E2 has non-zero occupation. For higher values of 
F, i.e., when the bands split, we determine the chemical potential fi = fi{n, F) corresponding 
to the fixed carrier density. The values of fi provided below represent zero field values. 
Shown in Fig. |4] is the dependence of the solution of the gap equation on the Zeeman 
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FIG. 4. (Color online) The solution of the gap equation plotted against the Zeeman splitting F 
for three different values of for a fixed value of the spin-orbit coupling constant a = O.leV-A. 
The value of a is large enough so that A remains continuous with F. For larger values of /u (black 
and green curves), A becomes inappreciably small for F > /i. For the red curve, however, A is 
appreciable ( 0.02 meV) for F > /x = 0.25 meV. Since to the right of this F it is possible to satisfy 
F^ > yU^+A^, the system is in a topologically non-trivial phase in this region. Therefore, somewhere 
above F = 0.25 meV (shown with an arrow), there is a topological quantum phase transition from a 
regular superconducting phase (to the left of the arrow) to a topologically non-trivial non-Abelian 
phase (to the right of the arrow). 

splitting for three different values of the chemical potential (i.e., three carrier densities) 
/i = 2.0meV, ImeV, 0.25meV. The Debye frequency is taken as ud = 25meV and the Rashba 
coupling is a = 0.1e\^-A, i.e., the system is characterized by a strong spin-orbit coupling. 
For these parameters the system is above the critical value of a for the discontinuous fall of 
A and hence A is now a continuous function of F. Before analyzing the plots in Fig. |4| let 
us remind ourselves that in order to satisfy the conditions for the non-Abelian s-wave phase 
(F^ > yU^ + A^), we need an appreciable A when F has crossed ~ jj. Coming back to Fig. |4]we 
note that, similar to Fig. |3] (black curve), A falls with increasing values of F. For the black 
and the green curves (higher /i) A becomes inappreciably small (though it is never zero) 
by the time F becomes ~ ^. However, for the red curve (/i = 0.25meV) there is a residual 
superconducting pair potential A ^ 0.02meV for F > /i, i.e., the system is in a topologically 
nontrivial non-Abelian phase. Moreover, as shown in Fig. [5} the magnitude of this residual 
s-wave pair potential increases with a, and thus can be increased by increasing the value of 
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FIG. 5. (Color online) Dependence of the superconducting pair potential A in the topological 
phase on the strength of Rashba spin-orbit coupling a. The system is characterized by a chemical 
potential at zero magnetic field [iq = 0.25meV (see the red curve in Fig. |4]for the dependence of 
A on r). The Zeeman splitting is fixed at F = O.SmeV. With these parameters the system is in a 
non-Abelian phase for a = O.leV-A. The superconducting pair potential increases monotonically 
with a. For a > 303mey-A(black arrow) the condition F^ > /u^ + A^ is no longer satisfied and the 
system undergoes another topological quantum phase transition to an ordinary superconducting 



phase (at large a). The inset shows the two excitation energies, Eq. (26), of which the smaller one 
(blue curve) vanishes at /c = at this a-tuned TQPT. 

the spin-orbit coupling. Therefore, for these parameter values, there is topological quantum 
phase transition (TQPT) when F crosses the critical value F^ = + (shown with an 
arrow in Fig. |4]). The TQPT separates a regular (non-topological) superconducting phase 
(F < Fc) from a topological non-Abelian superconducting phase (F > Fc). From our self- 
consistent mean field theory we find this TQPT to be continuous, that is, there is no change 
in A at the critical value of the Zeeman splitting. 



VIII. TQPT IN THE PROXIMITY INDUCED CASE 

An alternative and perhaps more robust way to create a topologically non-trivial non- 
Abelian superconductor is to induce a superconducting pair potential in a spin-orbit coupled 
semiconductor by proximity effeclP^. Ideally, for the proximity effect induced superconduc- 
tivity the pairing interaction resides in a parent s-wave superconductor such as Al or Nb 
while the quasiparticles of interest are confined to a two- or one-dimensional semiconductor 
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FIG. 6. (Color online) Quasiparticle gap Eg versus Zeeman coupling T for various values of spin- 
orbit interaction a. The strength of the spin-orbit coupling in the inset is such that a = 0.3 corre- 
sponds to 0.1 eV-A. The proximity induced pair potential and chemical potential are taken to be 
Ag// = O.bmeV and ^ = 0.0. The quasiparticle gap vanishes at the critical value Tc = y^^A^j^H-"/^- 
Above the critical point spin-orbit coupling opens a quasiparticle gap that is proportional to a in 
the small a limit. 

layer on the surface of the superconductor. The proximity effect has been shown to create 
a topological superconductor similar to the ones discussed above on the surface state of a 
topological insulatoi'^'^^ESl and also in a 2D semiconductor layei'^. 

Physically, the proximity effect arises from multiple Andreev reflections of electrons in a 
semiconductor that is connected to a superconductor by tunneling. For most realistic cases, 
there is no pairing interaction in the semiconductor. Thus, strictly speaking the supercon- 
ducting pair potential vanishes in the semiconductor and at first glance it appears that there 
is no superconductivity induced in the semiconductor. However, the superconducting order 
parameter defined by is found to remain non-zero in the semiconductor layer. 

Furthermore, the multiple Andreev reflections open a gap in the spectrum of quasiparticles 
that are localized in the semiconductor layer. The spectra of such quasiparticles can be 
shown to be identical to quasiparticles with an effective pairing potential in the semicon- 
ductor layei^^. Therefore from the point of the quasiparticle spectrum, which is the only 
property that is relevant to the definition of a topological superconductor, the proximity to 
a superconductor induces a superconducting quasiparticle gap in the semiconductor. 

The proximity effect can be induced by even weak tunneling between the semiconductor 
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and the superconductor. Therefore, the quasiparticle spectrum in the semiconductor does 
not affect the pairing potential in the superconductor significantly. Specifically, for the 
proximity induced superconductivity case, the self- consistency effects that were important in 
the discussions in the previous sections become insignificant. Furthermore, if the Zeeman 
potential is also induced by proximity effect from a magnetic insulator on the other surface of 
the semiconductor, there is no direct tunneling between the superconductor and the magnetic 
insulator and therefore no suppression of the order parameter in the superconductoi'^. Thus, 
in contrast to the discussions in the previous sections, where the Zeeman potential induced 
topological phase transition was accompanied by significant changes in the pair potential 
A, the pair potential A in the proximity induced case remains unaffected by the Zeeman 
splitting. 

The TQPT in both cases (the proximity induced case and the case when the pair po- 
tential is due to an intrinsic pairing interaction) can be characterized by the closing of the 
superconducting quasiparticle gap (shown in Fig. |6]) as the Zeeman potential is raised from 
F = past the critical value Fc = a/A^ + yU^. In the proximity induced case, A is the 
proximity induced effective pair potential and fi is the fermi energy in the semiconductor. 
The quasiparticle gap Eg{k) (minimum of Ei^2){k) in Eq. (26)) closes at = exactly when 
F passes through F^, (Fig. |6]), indicating the existence of a QPT even though the supercon- 
ducting pair potential A remains perfectly continuous. The quasiparticle gap for F > Fc 
shows a linear dependence on the spin-orbit coupling strength a at small cP. Here it is 
appropriate to mention a caveat for the case where the Zeeman potential is not proximity 
induced and instead induced by a magnetic fiel d^° * ^^ *^. In this case, the Zeeman potential 
also suppresses the superconducting pair potential in the parent s-wave superconductor. 
However, this effect can be small provided the (/-factor in the semiconductor is much larger 
than in the superconductor as is often the case in 2D electron systems. 



IX. DISCUSSION 

Topologically non-trivial non-Abelian superconductivity can be realized in two different 
classes of systems. In class (a), superconductivity is proximity induced on a semiconductor 
(in the form of a film or a wire) which has a strong spin-orbit coupling. In class (b), 
superconductivity arises from intrinsic attractive pairing interaction in a system which also 
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has a sizable spin-orbit coupling. In both cases a firm requirement for a phase transition 
from an ordinary superconducting phase to a topologically non-trivial superconducting phase 
is an externally imposed Zeeman splitting. The Zeeman splitting creates a gap in the spin 
orbit bands (Fig. [T]). When this gap is large (Zeeman splitting is comparable to the chemical 
potential) so the Fermi surface lies in only the lower band, it triggers a QPT at which the 
system goes from a regular superconducting phase (small Zeeman splitting) to a topolgical 
superconducting phase (large Zeeman splitting). This value of Zeeman splitting far exceeds 
the value at which an ordinary s-wave superconductor is known to lose superconductivity 
due to its inability to form spin-singlet zero-momentum Cooper pairs. As we have shown 
above, this is where the requirement of a sizable spin-orbit coupling is important to stabilize 
a topological superconducting phase. Below we recapitulate and discuss the main results 
first for the case where the superconductivity arises from an intrinsic pairing interaction 
followed by the much simpler case of superconductivity arising from proximity effect. 

To discuss the various phases and the quantum phase transitions we have divided the 
parameter space in two distinct regimes by the relative magnitude of fi with respect to all 
other energy scales in the problem. In the large /z regime the underlying system always has 
two Fermi surfaces irrespective of the magnitude of the Zeeman splitting F. In the absence 
of the spin-orbit coupling a, with increasing values of F it becomes increasingly difficult for 
the system to create spin-singlet s-wave pair potential at zero net momentum. Ignoring the 
possibility of Cooper pairs with non-zero net momentum, we find that when F crosses a 
critical value ~ A the system becomes non-superconducting at a first order QPT. At this 
transition the pair potential drops discontinuously to zero. By including a non-zero a we 



find that, surprisingly, there is always a non-zero solution of the gap equation Eq. (28). 
This is because with a ^ the individual bands can no longer be viewed as carrying a 
single spin component. Rather, both bands now carry a minority spin amplitude along 
with the majority component, which allows s-wave superconducting pairing even for large 
values of F. If a is small, 2akF < A(0) where A(0) is the value of the order parameter for 
zero Zeeman splitting, there is still a precipitous drop of A at Zeeman splitting F ~ A(0) 
(Fig. [s]). However, this is not a QPT, since, as already mentioned, A is never strictly zero 
in the presence of a non-zero a. When a itself crosses a threshold value, 2akF > A(0), the 
first-order-like drop of A as a function of F turns into a slower continuous decay (black curve 



in Fig. [3j). For high values of F > 2akF, A again decays exponentially with F. However, 
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this high field scale, comparable to the spin-orbit strength at the Fermi surface, may not be 
attainable in real systems. Consequently, the s-wave superconducting gap may never vanish 
with a Zeeman coupling in the presence of strong spin-orbit coupling. 

The regime of small fi is particularly important because of the possibility of a topological 
phase transition. In this case, the behavior of A with F for a = (first order QPT) and small 
a (precipitous drop of A with F) remain unchanged from the case with large /i. For small 
H, however, A(0) itself is small. Consequently a is always in the large spin-orbit coupling 
regime, 2akp > A(0). Therefore, for realistic values of a, A falls only gradually with F and 
strictly speaking is never zero (Fig. |4]). Let us now recall that for a TQPT from a regular 
s-wave superconductor to a topologically non-trivial superconductor the parameters need 



to satisfy the condition F > Fc = a/A^ -|- /i^. This implies that, for a robust non-Abelian 
phase, we require an appreciable A when F becomes > fi. From the red curve in Fig. |4| 
we note that for F ~ /i, A is still appreciable, A ~ .02meV, and thus a stable non-Abelian 
phase is in principle allowed. Moreover, as shown in Fig. [sj the value of A for large F (i.e., 
A in the non-Abelian phase) is directly related to the spin-orbit strength a and increases 
appreciably if a can be increased (as in a cold fermion system). Conversely, there is no 
non-Abelian phase (A = 0) if the system has no spin-orbit coupling. 

From our self-consistent mean field theory, we find the TQPT at F = F^ = a/AM-T^ 
to be continuous. By this we mean that the magnitude of A is continuous across this 
transition. At F = Fc the underlying system shifts from having two Fermi surfaces (F < Fc) 
to just one in the lower band (F > Fc). As shown in Ref. [Ij, for F > Fc a defect in the 
superconducting order parameter (e.g., vortex, sample edge) traps a unique zero energy 
bound state Majorana mode. Such a non-degenerate bound state solution is absent for 
F < Fc. The emergence of the topological bound state Majorana solution for F > Fc 
makes the transition a topological one. The exact location of the topological transition is 



indicated by the quasiparticle excitation energy Eg{k) (minimum of En2){k) in Eq. (26)) 
passing through zero. This happens at A; = exactly when F passes through Fc (Fig. |6]), 
indicating the existence of a QPT even though the superconducting order parameter A 
remains perfectly continuous. 

When s-wave superconductivity is proximity induced on a semiconductor, there is no 
self-consistent gap equation to be satisfied in the semiconductor. Thus there is no self- 
consistency effects that suppress the pair potential with the Zeeman splitting as discussed 
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above. In this case, the semiconductor simply 'inherits' the superconducting pair poten- 
tial and its quasiparticle spectrum is modified accordingly. For weak tunneling between 
the semiconductor and the superconductor layers, the quasiparticles in the semiconductor 
cannot significantly infiuence the pair potential in the host superconductor. Therefore, the 
self-consistency requirement as in the discussions above can be neglected. If the Zeeman 
potential is also induced by the proximity effect of a magnetic insulator from the opposite 
side of the semiconductor, there will be minimal effect of the magnetic insulator on the s- 
wave superconductor. If the Zeeman potential is induced by a parallel magnetic field, then 
the effect on the host superconductor will again be minimal provided the (^-factor in the 
semiconductor is larger than that in the superconductor. 



X. CONCLUSION 

To conclude, we have considered spin-orbit coupled systems with superconductivity aris- 
ing from either intrinsic on-site s-wave pairing interactions or from the proximity effect of 
an adjacent superconductor. In both cases, using BdG analysis of a postulated mean field 
Hamiltonian with an s-wave pair potential A, it has been showrP^that when an externally 
imposed Zeeman splitting crosses a critical value, there is a Majorana fermion mode at a 
vortex core. The required Zeeman splitting, F > Fc = y^/i^ + A^, seems to far exceed the 
value (F ~ A) above which an s-wave pair potential A is known to vanish. This gives rise 
to the conceptual question if the postulated pair potential in Ref. [1] and all the subsequent 
works on this system is spurious for F > Fc. If true this will indicate that the BdG result 
of the Majorana fermion at a vortex core for F > Fc, based on the postulated mean field 
is a spurious mathematical result with no physical connection. In this paper we have 
resolved this question by showing that in the presence of spin-orbit coupling the s-wave pair 
potential can never be made strictly zero by the application of a Zeeman potential. This is in 
agreement with a similar result derived previously in a different contextP. When the s-wave 
pair potential arises from an intrinsic local pairing interaction, our self-consistent analysis of 
the gap equation reveals that the decay of the pair potential with Zeeman splitting is more 
gradual in the presence of spin-orbit coupling, although for large enough Zeeman splitting 
the decay is exponential. Thus there can be a small but finite region in the parameter space 
(RHS of arrow, which indicates a TQPT, on Fig. 4) where a topologically non-trivial super- 
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conducting phase can be realized. When the s-wave pair potential is proximity induced on 
a semiconductor, there is no requirement of satisfying the self-consistent gap equations. In 
this case, the pair potential is simply 'inherited' from the adjacent superconductor. Thus in 
this case the topologically non-trivial phase is much more robust than the other case where 
it is due to intrinsic pairing interactions. 

Two final comments are in order here. For long- but finite-range pairing interactions (as 
opposed to local interactions as in this paper) it is well known that the spin-orbit interaction 
mixes s-wave and p-wave pair potentialJ^. In this case it may appear that superconductivity 
can evade the CC limit merely because the p-wave part of the pair potential can survive 
the strong Zeeman field, even if the s-wave part cannot. It is, however, incorrect to ascribe 
the existence of the topological superconducting state at large T to this effect. As we have 
shown in detail in Ref . [1] , the topological state owes its existence solely to the survival of the 
s-wave part of the pair potential. (The Pfaffian topological invariant discussed in Ref. [1] is 
completely insensitive to the p-wave part.) It is precisely to isolate and eliminate the effect 
of the mixed p-wave pair potential that in this paper we confined ourselves to a strictly local 
pairing interaction. The existence of the topological state at high Zeeman fields is strictly 
due to the survival of the s-wave pair potential, the physics of which is discussed in this 
paper and also summarized in the concluding paragraph of the introduction. 

But this work is not just an academic resolution of the question of the survival of an s- 
wave pair potential in the presence of a strong Zeeman field. It also proves that all properties 
of the topological state when superconductivity is proximity induced continue to hold even 
when superconductivity is due to local s-wave pairing interactions. This result is directly 
relevant to the case of an s-wave Feshbach cold atom system. Note that in this case A cannot 
just be assumed in the BdG equations (as in the case of the proximity effect frameworlP^, 
but has to be calculated from the gap equations as in the present paper. 
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